We investigate strategies to numerically integrate closed lines and surfaces that are implicitly defined by level sets (iso-contours) of continuously differentiable toroidally symmetric functions. The “grid-transform” approach transforms quantities given in non-surface-aligned coordinates onto a numerically constructed surface-aligned grid. Here, line and surface integrals, as well as so-called ...